Bond disorder and breakdown of ballistic heat transport in the spin-1/2 
antiferromagnetic Heisenberg chain as seen in Ca-doped SrCu02 
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We study the impact of a weak bond disorder on the spinon heat transport in the 5 = 1/2 anti- 
ferromagnetic (AFM) Heisenberg chain material Sri_xCaxCu02. We observe a drastic suppression 
in the magnetic heat conductivity Kmag even at tiny disorder levels (i.e., Ca-doping levels), in stark 
contrast to previous findings for Kmag of S = 1/2 two-dimensional square lattice and two- leg spin- 
ladder systems, where a similar bond disorder has no effect on Kmag. Hence, our results underpin 
the exceptional role of integrability of the S = 1/2 AFM Heisenberg chain model and suggest that 
the bond disorder effectively destroys the ballistic nature of its heat transport. We further show 
that the suppression of ftmag is captured by an effective spinon-impurity scattering length, which 
exhibits the same doping dependence as the long-distance exponential decay length of the spin-spin 
correlation as determined by density-matrix renormalization group calculations. 

PACS numbers: 75.40.Gb, 66.70.-f, 75.10.Kt, 75.10.Pq 



I. INTRODUCTION 

The transport properties of the integrable one- 
dimensional (ID) 5 = 1/2 antiferromagnetic (AFM) 
XXZ chain model are attracting considerable atten- 
tion because anomalous spin and heat transports have 
been predicted.—^— Rigorous predictions concern the heat 
transport of the model in the isotropic (Heisenberg) case. 
It is known to be ballistic as the consequence of integra- 
bility and fundamental conservation lawsj^"— This means, 
a divergent heat conductivity is expected. However, in 
experimental realizations of this model, the observable 
spinon heat conductivity K,nag always has to be finite 
since extrinsic scattering processes due to defects and 
phonons are inherent to all materials and mask the in- 
trinsic behavior of a spin chain. Nevertheless, a very 
large Kmag bas been observed in a number of cuprates 
that reahze 5 = 1/2 spin chains>i2rJ^ Especially, the ma- 
terial SrCu02 is considered an excellent realization of the 
ID- AFM 5 = 1/2 Heisenberg model (HM)J^-^ For this 
material, we recently reported quasi ballistic spinon heat 
transport with mean-free paths > 1 pm for samples of 
extraordinary purityJ^ 

Recently, several theoretical studies focused on the 
consequences of disorder for transport in ID systems 
with controversial resultsi^"— Experimental studies are 
scarce and, so far, only concern heat transport studies on 
the ID- AFM 5 = 1/2 HM with site disorder (diagonal 
disorder), which cuts the chain into finite segments and, 
thus, rather trivially leads accordingly to a suppression of 
transportiii Here, we study the more subtle bond disorder 
(off-diagonal disorder) , which leaves the chains in the ma- 
terial intact and which recently has been investigated in 
nuclear magnetic resonance (NMR) experiments i^iS^ We 



induce this subtle type of disorder in SrCu02 by system- 
atically substituting isovalent Ca for Sr in tiny amounts. 
This off-chain lattice doping by small Ca^+ ions induces 
a local lattice distortion which must result in a mod- 
ulation, i.e., disorder, of the magnetic exchange con- 
stant J. In contrast to other doping schemes where the 
dopants occupy Cu sites within the magnetic structures 
and, thus, generate site disorder fiii^i^ for this bond dis- 
order, one generally expects only subtle changes in the 
magnetic transport properties of the system. In fact, pre- 
vious studies show that the magnetic heat transport of 
the 5=1/2 AFM square lattice^^ and the two-leg spin 
ladder systems?^ whose underlying spin models are non 
integrable, is unaffected by such kinds of disorder. Radi- 
cally different from these findings, in Ca-doped SrCu02, 
we observe a severe suppression of Kmag already at tiny 
doping levels (~ 1%). Hence, these data suggest that the 
disorder-induced departure from integrability efficiently 
destroys the ballistic nature of heat transport in the ID- 
AFM 5=1/2 HM. In the framework of a kinetic model, 
we show that this suppression is captured by an effec- 
tive spinon-impurity term in the spinon mean-free path, 
which surprisingly exhibits the same doping dependence 
as the long-distance decay length of the spin-spin corre- 
lation as determined by density-matrix renormalization 
group (DMRG) calculations. 



II. EXPERIMENTAL DETAILS 

We have grown centimeter sized single crystals of 
Sri_^Ca^Cu02 with a; = 0, 0.0125, 0.025, 0.05, 0.1 using 
the traveling-solvent floating zone method<^ As starting 
materials we used CuO (99.99% purity), SrCOs (99% 
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purity) and CaCOs (99% purity). Additionally, crys- 
tals with x=0, 0.0125 were grown with all starting pow- 
ders of 4N (99.99%) purity. The crystallinity and the 
doping profile were checked under polarized light and 
by energy dispersive x-ray spectroscopy. A structural 
refinement was performed for samples of SrCu02 and 
Sro.9Cao.iCu02 in a single-crystal diffractometer. For 
the doped material, this yielded a reduction in the cell 
volume by 1% with respect to the undoped compound 
in agreement with literature values.— Furthermore, we 
found a change in the Cu-O-Cu bond angle and the 
bond distance of roughly 0.3%^. Since these data rep- 
resent averages over the entire crystal volume, we esti- 
mated the local Ca-induced variation in the lattice from 
density-functional calculations using the code Quantum 
EspressOj24 a 2 X 2 X 2 superceU of SrCu02 doped with 
10% Ca was relaxed with minimal symmetry assump- 
tions. This yielded a Cu-O-Cu bond length variation that 
is approximately the same as the doping induced change 
in the mean Cu-O-Cu distances measured with the single- 
crystal diffractometer. This consistency corroborates the 
presence of a bond disorder in the material. For the 
transport measurements, rectangular samples with typi- 
cal dimensions of (3 ■ 0.5 ■ 0.5) mm'^ for measurements of 
the heat conductivity along the principal axes (k^, k^, 
Kc) were cut from the crystals for each doping level with 
an abrasive slurry wire saw. The longest dimension was 
taken parallel to the measurement axis. Measurements 
of the thermal conductivity as a function of temperature 
T in the range of 7 — 300 K were performed with a stan- 
dard four probe technique^ using a differential Au/Fe- 
Chromel thermocouple for determining the temperature 
gradient. 



III. EXPERIMENTAL DATA 

Figure [1] shows our experimental results for Kc and 
Ka of Sri_2;Caa;Cu02 for all doping levels. In the pris- 
tine material {x = 0), the heat conductivity parallel 
to the spin chains Kc is strongly enhanced with respect 
to the purely phononic heat conductivity perpendicular 
to the chains Ka- This is a result from quasi ballis- 
tic spinon heat transport in the chainsJ^ More specif- 
ically, the spinon heat conductivity Kmag, which domi- 
nates the total (spinon plus phonon) heat conductivity 
along the c axis leads to a broad peak at low temper- 
atures (kc « 830Wm"iK"i at T « 28 K) followed by 
a strong decrease toward a still significantly enhanced 
value (^ 50Wm~^K^^) at room temperature. 

Upon doping the material with Ca, the thermal con- 
ductivity perpendicular to the spin chains (kq) is gradu- 
ally suppressed. This is the typical expectation of in- 
creased scattering by phonons off defectsi^ Since the 
Ca^"*" impurities possess a smaller ionic radius and 
mass as compared to Sr^+, they act as defects for the 
phonons4S A much more dramatic change upon doping 
is observed in Kc- Already at the lowest doping level of 
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Figure 1: (Color online) Thermal conductivity parallel to the 
spin chain (kc) for Sri_a;Caa;Cu02 at different doping levels. 
(Inset) Thermal conductivity for the same doping levels per- 
pendicular to the spin chain (na). 



1.25% Ca, Kc is suppressed strongly as compared to that 
of the pristine material, and the overall curve shape of 
Kc{T) is changed completely. The peak at low tempera- 
tures is now significantly smaller (k^ « 320Wm~^K~^), 
is much sharper, and is shifted toward lower tempera- 
tures (T 18 K). At about 50 K, Kc shows a kink, above 
which, it decreases much weaker and approaches Kc of the 
undoped compound for even higher temperatures. Upon 
increasing the Ca doping, the absolute value of the peak 
decreases, although its position stays constant. The kink 
is gradually shifted to higher temperatures, and all of the 
curves approach the curve of the undoped compound at 
room temperature. This apparent saturation behavior is 
most visible at 10% Ca doping, where a small maximum 
is followed by a practically constant thermal conductiv- 
ity. Nevertheless, the anisotropy between Ka and Kc is 
still considerable since at 300 K, Ka « 4Wm~-'^K~^ and 
Kc R:! 40Wm-iK-i4l 



IV. ANALYSIS 

The strong suppression of Kc implies that the Ca im- 
purities lead to a very strong suppression of the spinon 
heat conductivity Kmag already at very low doping levels 
{x « 0.01) because the doping scheme apparently only 
causes a moderate suppression of the phononic heat con- 
ductivity, as seen in k^. Qualitatively, this suggests that 
the doping-induced bond disorder leads to a significant 
deviation in the spin system in Srx_a;Ca2;Cii02 from tliG 
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Figure 2: (Color online) Kmag for the different doping levels 
as a function of temperature. (Inset) Enlarged portion of the 
main plot shows the different doping levels in more detail. 
The black dashed line indicates the shift in the maximum to 
higher temperatures for higher dopings. For 10% Ca doping, 
an extrapolation is shown without a spin gap (solid yellow 
line) and with a spin gap of 50 K (dotted orange line). The 
shaded areas show the uncertainty of the estimation of Kmag 
due to the phononic background. All curves are shown only 
down to a temperature for which the uncertainty of the esti- 
mation of Kmag is Still reasonably small. 

S ~ 1/2 Heisenberg chain model, and, thus effectively 
destroys the ballistic nature of the heat transport in the 
material. One might conjecture that the connected de- 
parture from integrability can be captured by an effective 
scattering process that describes the observed suppres- 
sion of K„iag in our heat transport experiments. In order 
to investigate this notion further, we analyze our data 
by extracting Kmag and calculating the spinon mean-free 

path Zmag- 

A. Spinon heat conductivity 

The first step of such an analysi o^^i^^ consists in esti- 
mating the phononic part of Kc via Kc,ph ~ ^a- Then, the 
spinon thermal conductivity is given by K„iag = Kc — k,/^. 
The doping-dependent evolution of Kmag is shown in 
Fig. [2] At low temperatures, Kmag of the undoped com- 
pound shows a large maximum (Kmag ~ 665Wm~^K~^ 
at T « 36 K) , which is followed by a steep decrease upon 
approaching room temperature where the values of the 
curve nearly become constant (K^ag ~ 42Wm~^K~^). 
Doping of 1.25% Ca leads to a severe suppression of 
Kmag at low temperatures with a broadening of the max- 
imum (^ 96 Wm~-'^K~^) and a shift to 69 K. Up to room 
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Figure 3: (Color online) Imag of Sri_3;Ca3;Cu02 for different 
levels of doping. Toward low temperatures, the error in the 
estimation of Kmag increases due to an increase in Kph. Thus, 
the values of Zmag are shown only at temperatures above which 
the error in Imag is reasonably small. The solid lines were 
calculated according to Eq. ([2]). 
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Table I: Fit parameters for the mean-free paths according 
to Eq. ([2|. Additionally, Iq is given in units of the lattice 
constant c. 



temperature, the curve approaches that of the undoped 
compound. Further increasing the doping continues to 
decrease Kmag at low temperatures and shifts the increas- 
ingly broadened maximum to even higher temperatures. 
For temperatures T > 200 K, the values of Kmag approach 
those of the undoped compound but remain smaller. 

B. Spinon mean free path 

For each doping level, we use the derived Kmag data 
to calculate the spinon mean-free path /mag (see Fig. [3]) 
according toi^ii^iiiii^ 

with Ng ~ A/ab as the number of spin chains per unit 
area. For all doping levels, /mag (7^) decreases strongly 
with increasing temperature which is the signature of 
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spinon-phonon scattering J^iis However, upon increasing 
the doping level, Zmag exhibits a systematic reduction, 
and the overall temperature-dependent change at a given 
doping level becomes less pronounced, suggesting an in- 
creased importance of spinon-defect scattering We test 
this notion by applying Matthiesen's rule and decompos- 
ing ^mag into two respective contributions, i.e., ^magC^) = 
^o"^ + ^sp{T)~'^. Here, Iq accounts for temperature- 
independent spinon-defect scattering, whereas, lsp(T) 
corresponds to T-dependent spinon-phonon scattering. 
By using l~p oc T cxp{~T*/T) with a characteristic 
phonon energy scale T* on the order of the Debye 
temperature,— liii^ we have 

l7nisiT) = lo' + A.Texp {^T:/T) , (2) 

with a proportionality constant Ag . Equation ([2]) is used 
to fit the ^mag data in Fig. [31 where T* first was de- 
termined for the pure compound and then was used as 
a constant when fitting the data for finite-doping lev- 
els. As can be seen in Fig. [31 the fits describe the data 
very well (see Table [J for the fit parameters)^. A re- 
markable result is that As also is constant nearly for 
all doping levels with a root-mean-square deviation of 
around 10% of the arithmetic mean value. Thus, the 
effect of the bond disorder induced by the Ca doping 
is described well by only one doping-dependent parame- 
ter that describes the spinon-defect scattering Iq and one 
temperature-dependent spinon-phonon scattering mech- 
anism, which is independent of doping. 

Alternatively, for the afore discussed scattering, one 
might try to attribute the suppressed K^ag to a disorder- 
induced depletion of thermally excited quasiparticles. In 
fact, recent NMR measurements provided evidence for a 
spin gap A/fcs w 50 K for 10% Ca dopingi^ However, 
the large J makes the effect of this gap on the transport 
tiny. To illustrate the negligible influence, the inset in 
Fig. [21 shows an estimate of the impact of a spin gap^ 
on Kinag at 10% Ca doping. 

Having established that the suppression of Kmag, in- 
duced by the bond disorder, indeed can be described by 
an effective spinon-defect scattering process, in Fig. [H 
we investigate how the corresponding scattering length 
Iq develops as a function of doping. As can be inferred 
from the figure, Iq decreases systematically with doping, 
following ^0 = ai/a;-|-/iim, with the scattering strength ai 
and the offset /lim. Note that, while Iq ^ 1/x is consistent 
with previous findings for spin-defect doped chains, two- 
leg spin ladders and 2D square latticesjiii^i^S the finding 
of a finite offset ln^ is rather unexpected. Qualitatively, 
such a finite offset means that the effect of the bond dis- 
order, in terms of defect like scattering, is strong at very 
small degrees of disorder but becomes increasingly unim- 
portant at larger degrees of disorder. Hence, the effect 
of the bond disorder is significantly different from that 
of strong site defects (diagonal disorder), which cut the 
chains into finite segments 
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Figure 4: (Color online) lo and ^ as a function of doping 
X (bottom axis) and disorder parameter D (top axis), re- 
spectively. The quantities are fitted by Zo = ai/x + 
(ai « 4.3, hm. ~ 29) and ^ = aa/D-fCum ("2 ~ 16, « 26). 
It is worth noting that only a simple scaling factor for D is 
necessary to reach an almost perfect agreement between both 
curves. 



C. Spin-spin correlation 

In order to obtain deeper insight into the effect of bond 
disorder on the spin dynamics of the system, we investi- 
gated its influence on the spin-spin correlation as a func- 
tion of the distance z. The model Hamiltonian is written 
as 

H = J2iJ + De^)S,■S,+^, (3) 

i 

where Si is a spin-^ operator at site i, J is the nearest- 
neighbor exchange interaction without bond disorder, 
Ei is defined by a box probability distribution V{ei) = 
0(1/2 — \ei\) with the step function 0{x), and the dis- 
order strength is controlled by D. The spin-spin cor- 
relation functions (^Si ■ Si+z^ are investigated using the 

DMRGi^ We study chains with 1024 sites keeping 
m = 2000 density-matrix eigenstates. Here, the open- 
boundary conditions are used, and the distance z is cen- 
tered at the middle of the system. We calculate the corre- 
lation functions by randomly sampling 300 realizations of 
Visi) = 6(1/2— \ei\) and taking an average of the results 
for each D and z. For D = 0, the spin-spin correlation 
follows a 1/z law, which is expected in the ballistic case. 
For D > the long-distance part can be fitted just by 
exp(— z/^), with the free parameter ^. Since the com- 
mutator of the heat current operator and the spin chain 
Hamiltonian is no longer zero for D > 0, ballistic heat 
transport no longer is expected. If one interprets the 
deviation in the spin-spin correlation from the algebraic 
decay (l/z) as a measure for the departure of the system 
from the clean one (i.e. from integrability) , the argument 
of the fitted exponential may fix the crossover, with z < ^ 
being ballistic and z > ^ being diffusive. Plotting ^ as a 
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function of D (cf. Fig.[3]) strikingly reveals that ^ depends 
in the same manner on D as does on x. Note, in par- 
ticular, that the calculated offset quantitatively almost 
perfectly matches that of the experiment, i.e., ^ « in 
the strong disorder or large doping limit. Furthermore, 
one may apparently assume D oc x. Hence, in the frame 
of this simple model, the long-distance decay length ^ of 
the spin-spin correlation can be interpreted as a limit for 
the effective spinon-defect scattering length Iq. 

V. CONCLUSION 

In conclusion, the doping-induced bond disorder in 
Sri_a;Ca2;Cu02 causes a severe suppression of the spinon 
heat conductivity which is consistent with a disorder- 
induced departure of the underlying spin model from 
integrability and, thus, the destruction of ballistic heat 
transport. This interpretation is corroborated by the fact 
that the observed suppression apparently is a character- 
istic feature of this ID-AFM S = 1/2 HM because no 



effect of bond disorder has been observed for a non in- 
tegrable S = 1/2 AFM square lattice and two-leg ladder 
systems. Using a simple kinetic model, we have shown 
that the suppression of Kmag is described well by an effec- 
tive spinon-defect scattering length, which can be related 
to the long-distance decay length of the spin-spin corre- 
lation as calculated by the DMRG method. 
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We find a bond length and angle of 3.910 A (3.896 A) and 

175.5° (174.9°), respectively, for a; = (x = 0.1). 

We also checked the thermal conductivity perpendicular to 

the chains along Kb. This gives very similar results apart 

from a small anisotropy, which is already present in the un- 

doped compound and does not change in magnitude upon 

doping. 

Note, that measurements of Ate with 2N and 4N purity of 
1.25% Ca doping did not show any significant difference. 
Therefore the results for 2.5%, 5% and 10% Ca doping were 
obtained using crystals of 2N purity. 

Considering the small differences between Ka and Kt of 
~ 15% compared to the overall anisotropy of more than a 
factor of 5 between Kc and K,a, the validity of this method 
seems to be justified. 

Modeling the spinon scattering by thermally excited opti- 
cal phononsi^i^ achieves a fit of similar quality and com- 
parable Iq but with different T* . 



